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Abstract 

The KASCADE experiment measures extensive air showers induced by cosmic rays 
in the energy range around the so-cahed knee. The data of KASCADE have been 
used in a composition analysis showing the knee at 3-5 PeV to be caused by a steep- 
ening in the hght-element spectra [1]. Since the apphed unfolding analysis depends 
crucially on simulations of air showers, different high energy hadronic interaction 
models (QGSJet and SIBYLL) were used. The results have shown a strong depen- 
dence of the relative abundance of the individual mass groups on the underlying 
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model. In this update of the analysis we apply the unfolding method with a differ- 
ent low energy interaction model (FLUKA instead of GHEISHA) in the simulations. 
While the resulting individual mass group spectra do not change significantly, the 
overall description of the measured data improves by using the FLUKA model. In 
addition data in a larger range of zenith angle are analysed. The new results are 
completely consistent, i.e. there is no hint to any severe problem in applying the 
unfolding analysis method to KASCADE data. 



1 Introduction 



Due to the rapidly falling intensity with increasing energy, cosmic rays of ener- 
gies above 10^^ eV can be studied only indirectly by observations of extensive 
air showers (EAS) which are produced by the interactions of cosmic parti- 
cles with nuclei of the Earth's atmosphere. The observation of a change of 
the power law slope [2] of the size spectrum of EAS and consequently of the 
all-particle energy spectrum at ^ 3 • 10^^ eV 50 years ago has not yet been 
convincingly explained [3]. Several theories for the origin of the knee predict 
different knee positions for particles of different primary mass. Therefore the 
energy spectra of single elements or at least mass groups are of considerable 
interest. Indeed, recent analyses [4,1] find a steepening in the energy spectra 
of the light components in the knee region. Whereas the measurement method 
of detecting air showers alleviates statistical problems, one has to rely on the 
results of simulations and the description of hadronic interactions while re- 
constructing the properties of the primary particles. Since the required energy 
and important kinematic regions of these interactions are beyond the range of 
collider or fixed target experiments, the interaction models used are uncertain 
and differ in their predictions. On the other hand, a thorough analysis of EAS 
data offers the opportunity of testing [5] and improving the validity of these 
hadronic interaction models. 

The KASCADE-Grande experiment [6], located on site of the Forschungszen- 
trum Karlsruhe (Germany), is designed to measure EAS in the energy range 
between 0.5 PeV and 1 EeV. The installation consists of the original KAS- 
CADE [7] experiment and an extension by the Grande array, covering an 
effective area of 0.5 km^. 

The data analysis pursued by the KASCADE collaboration invokes an un- 
folding procedure of the two-dimensional shower size spectrum (total elec- 
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tron number vs. muon number) into energy spectra of five individual mass 
groups [1]. Despite the success of this method for the reconstruction of the 
shape of the spectral forms, a strong dependence of the result for the elemen- 
tal abundances on the interaction model underlying the analysis was found. 
Also, an insufficient description of the measured data by the employed simu- 
lations has been demonstrated. While in [1] the results for the analysis based 
on two different high energy interaction models were compared, we give in 
the present paper an update of the composition analysis concentrating on the 
influence of the low energy interaction model on the result. Furthermore, the 
analysis is repeated for measured data of different zenith angle intervals, thus 
testing the consistency of the analysis. 



2 Composition analysis of the KASCADE shower size spectrum 



Starting point of the analysis is the so-called two-dimensional shower size 
spectrum, i.e. the number of measured EAS depending on the electron number 
Ig and the truncated muon number Ig A^^^ (number of muons with shower 
core distances between 40 m and 200 m). In Fig. 1 this spectrum is given for 
showers inside the KASCADE array for three different ranges of inclination 
angle: EAS between 0° and 18^ between 18° and 25.9°, and between 25.9° 
and 32.3°, respectively. 

The content Nj of each histogram cell j can be written as 

'"T' rll 



C is a normalizing constant (time, aperture), and the sum is carried out over 
all primary particle masses A. The functions pa — Pa{^S Agj, Ig A^j | Ig E) give 
the probability for an EAS with primary energy E and mass A to be measured 
and reconstructed with shower sizes N^j and Nj[j. The probabilities pa include 
shower fluctuations, efliciencies, and reconstruction resolution. For reasons 
of clarity integration over solid angle and cell area is omitted in Eqn. 1. In 
case of KASCADE data pa is dominated by the shower fluctuations, whereas 
reconstruction systematics play an inferior role [1]. The data range was chosen 
in a way to minimize influences from inefficiencies. 

Adopting this notation the two-dimensional shower size spectrum is regarded 
as a set of coupled integral equations. In the analysis the primary particles 
H, He, C, Si, and Fe are chosen as representatives for five mass groups. The 
corresponding probabilities pa are determined by Monte Carlo simulations 
using CORSIKA [8] and a detailed GEANT [9] based simulation of the exper- 
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Fig. 1. Two-dimensional shower size spectra as measured by KASCADE for three 
different zenith angular ranges, where the first one shows similar data as used for 
the analysis described in [1]. The effective time of the measurements amounts to 
840 days. 

iment. To solve the equation system for the mass group energy spectra, the 
unfolding algorithm proposed by Gold [10] is applied. Details of the selection, 
reconstruction, and the analysis can be found in Ref. [1]. In particular, in that 
paper the determination and definition of the systematic uncertainties and 
inaccuracies of the analysis chain as well as the applicability of the chosen 
unfolding algorithm method is discussed in detail. 



3 Influence of low energy hadronic interaction model 

In the original analysis [1] the probabilities pa were determined using the high 
energy hadronic interaction models QGSJet [11] (2001 version) and SIBYLL [12] 
(version 2.1) in the simulations. For both cases the GHEISHA [13] code was 
used for interactions with laboratory energy of < 80 GeV. GHEISHA has 
been widely used as low energy hadronic model in air shower simulations in 
the last decade. 

Though the high energy interaction models are believed to have a larger effect 
on the KASCADE observables, the influence of the low energy hadronic model 
ought to be investigated. In particular, the number of secondary muons at sea 
level could be affected because these are decay products of low energy charged 
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Fig. 2. Distribution of electron (left) and truncated muon (right) shower size for iron 
induced showers of 1 PeV using GHEISHA and FLUKA as low energy interaction 
models and in both cases QGSJet as high energy interaction model. 




Fig. 3. Measured two-dimensional shower size spectrum ( — 18°) together with lines of 
the most probable values for proton and iron induced showers for both simulations. 

mesons. However, only at large distances from the core (where KASCADE is 
not sensitive), the contribution from low energy pions attains importance [17]. 
It is known [14,15] that the GHEISHA model suffers from deficiencies in han- 
dling the reaction kinematics properly leading to, e.g., a notably fiatter lat- 
eral distribution than the FLUKA model. In addition, GHEISHA produces 
too many pions as compared to accelerator data, while the FLUKA (version 
2002.4) [16] package is in line with the measured data. 



Therefore, in the more recent analysis of the KASGADE data GHEISHA is re- 
placed by the FLUKA model (preserving QGSJet version 2001 as high energy 
interaction model). Since FLUKA is known to describe the accelerator data to 
higher interaction energies, the change to the high energy interaction model 
is now at 200 GeV, i.e. FLUKA replaces in the intermediate energy range also 
the QGSJet model. 

When performing the simulations including detector response and the stan- 
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dard KASCADE reconstruction, the differences between these two models 
are found to be rather smaU, as can be seen in Fig. 2. As an example, 
the distributions of electron and truncated muon sizes for nearly vertical 
1 PeV iron induced showers are displayed, one simulated with the combina- 
tion QGSJet/GHEISHA, the other with QGSJet/FLUKA. In case of FLUKA 
simulations the shower size distribution is shifted by AlgA^g ^ 0.01 2% 
for both, proton and iron induced showers) towards larger electron numbers, 
and by AlgNj[ ^ 0.02 towards 5% for iron induced EAS, but only 2% for 
proton induced EAS) smaller muon numbers, i.e. the differences are larger for 
the muon number than for the electron number0 For illustration in Fig. 3 
the most probable values of the simulated electron-muon number correlation 
for the two models and for different primaries are displayed on top of the 
data distribution. The differences are much smaller than compared to those 
observed for different high-energy interaction models (see Fig. 22 in ref. [1]). 



10^ 



10 



O Proton, GHEISHA 
□ Helium, GHEISHA 
A Carbon, GHEISHA 



10^ 



• Proton, FLUKA 
■ Helium. FLUKA 
A Carbon, FLUKA 



> 

o 



★ Silicon, GHEISHA A Iron, GHEISHA 
^ Silicon, FLUKA A Iron, FLUKA 



10^ 



- 10^: 
^ 10: 



* 

^ 1 1 i ? 1 1 



10' 10' 
primary energy E [GeV] 



10' 



If 



10^ 10*^ 
primary energy E [GeV] 



Fig. 4. Comparison between QGSJet/FLUKA based results and QGSJet/GHEISHA 
based results for the energy spectra of H, He, and C (left) and of Si and Fe 
(right). Shaded bands correspond to estimates of the systematic uncertainties for 
the QGSJet/GHEISHA solutions. 



This behavior is expected in view of the above mentioned differences between 
GHEISHA and FLUKA: GHEISHA produces a higher pion multiphcity but 
flatter muon lateral distributions, or in other words, FLUKA predicts fewer 
muons in the range of 40 m - 200 m core distance. 



After carrying out the complete unfolding analysis, the results for the FLUKA 
case differ only little from the GHEISHA case, as might have been expected 
from the small differences in particle numbers. In Fig. 4, left panel, the results 



^ Investigations of these kinds by using SIBYLL instead of QGSJet as high energy 
interaction models lead for 1 PeV to a similar general behavior (^ 1% higher elec- 
tron number in case of FLUKA and a smaller muon number), but interestingly for 
SIBYLL there is a ^ 4% effect in muons for primary protons and only ^ 2% for 
primary iron. 
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Fig. 5. Distribution of individual x] ^he data range for the QGSJet solution. 
The low energy interaction model used is GHEISHA in the left panel (picture taken 
from [1]) and FLUKA in the right panel. 

for the energy spectra of H, He, and C using GHEISHA and FLUKA are com- 
pared with each other. It is a specific characteristic of the QGSJet model that 
the analysis of the KASCADE data always results in a Helium dominated 
composition. Differences between the two solution sets are small, especially 
when compared to the systematic uncertainties imposed by the unfolding pro- 
cedure [1] (represented by the shaded bands in the figure, which are displayed 
for the QGSJet/GHEISHA solution only, but which are of the same order for 
all the discussed solutions). This also holds for the corresponding spectra of 
Si and Fe, which are compared to each other in the right panel of Fig. 4. Here, 
the infiuence of replacing the low energy interaction model on the result is 
larger than for the light elements. The relative abundance of the heavy ele- 
mental groups increases slightly when using FLUKA for the unfolding. This is 
understandable taking into account the shift of the simulated electron-muon 
numbers of primary iron towards the bulk of the data distribution (Fig. 3). 
Finally, the results for the all-particle energy spectrum (which is the sum of 
the individual mass group spectra) are displayed in Fig. 9 together with re- 
sults discussed in the next section. No systematic difference in the all-particle 
spectrum could be observed by changing the low energy interaction model. 

To test the quality of the solution the obtained mass spectra are folded forward 
and compared to the measured data by a x^-test [1]. Noticeably, the Xdof 
parameter improves considerably from 2.38 to L 34 by going from GHEISHA to 
FLUKA. In Fig. 5 the individual x] distributions for both cases are displayed, 
suggesting that FLUKA describes the correlation of muon to electron number 
in air showers better than GHEISHA. 

To summarize, using the FLUKA model instead of the GHEISHA model seems 
to have no significant effect on the overall picture of the solution, but a slightly 
better description of the data is achieved. This improvement in describing the 
KASCADE data is not surprising considering the fact that FLUKA has been 
tuned to provide a good description of recent accelerator data. 
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4 Investigation of different KASCADE data sets 



In the analysis described so far, only EAS reaching the detector with zenith 
angles below 18° were included in the composition studies. The analysis of 
more inclined shower data can serve as consistency check. Due to the lim- 
ited reproduction of the measured observable correlations by the models, one 
cannot expect to obtain identical results for the energy spectra as compared 
to the vertical data set. Nevertheless, large differences between the solution 
sets for different zenith angle ranges could indicate a severe problem in ei- 
ther the simulation chain or the performed analysis technique. In particular, 
characteristics of the predictions of the attenuation of the electromagnetic and 
muonic shower components in the atmosphere can be tested. For this reason 
the analysis based on the QGSJet/FLUKA simulations were repeated for two 
additional data sets, the first one containing EAS with zenith angles ranging 
from 18° to 25.9°, the second one from 25.9° to 32.3°; thus covering the same 
acceptance on the sky. Data selection is identical to the one applied to the 
vertical data set, i.e. the same runs in the same period and the same cuts (ex- 
cept for reconstructed zenith angle). The two-dimensional measured shower 
size spectra are included in Fig. 1. 

Again, first the shower fluctuations are compared to get an impression on 
the differences in the data sets. In Fig. 6, as an example the distributions 
of electron and truncated muon sizes for 1 PeV proton induced showers are 
displayed for all three angular ranges. Going from nearly vertical to ^ 30° in- 
clined showers, the shower size distribution is shifted by A Ig Ag ?^ 0.5 towards 
smaller electron numbers (corresponding to an attenuation of the electron size 
relative to vertical EAS by ^ 70%), and by AlgNj[ 0.1, (attenuation of 
muons relative to vertical EAS by ^ 20%) towards smaller muon numbers. 
The widths of the distributions vary only slightly for the different angular 
ranges. Qualitatively, these differences are reflected in the two-dimensional 
shower size spectra shown in Fig. 1, considering the coverage of the measured 
distribution in particle numbers ill 

By applying the full unfolding analysis for the resulting mass group spectra 
only small differences can be detected (Fig. 7). The display of the systematic 
error bands is omitted in this figure for a better visibility. They are of the same 
order as those shown in Fig. 4 and therefore cover all three proton spectra, 
in particular at high energies. In the case of the proton spectra, systematic 
differences are present for energies above the proton knee. As can be seen, the 
knee in the proton spectrum gets less pronounced, i.e. the change of index 



^ The variation of the particle numbers with inclination is in first order independent 
of hadronic interaction models as it is due to attenuation in the atmosphere, which 
is handled by the CORSIKA code itself rather than by the interaction models. 
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Fig. 6. Distribution of electron (left) and truncated muon (right) shower size for 
proton induced showers of 1 PeV for different zenith angular ranges. The EAS 
are simulated with the QGSJet/FLUKA model combination. The lines show the 
parameterization of the distributions used in the analysis [1]. 
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Fig. 7. Results for the energy spectra of primary H and He, based on the analysis 
of EAS data originating from different zenith angle ranges (QGSJet/FLUKA). The 
display of the systematic uncertainties is omitted for reasons of clarity. 

decreases with increasing zenith angle. This may be related to the fact that 
shower fluctuations increase significantly for inclined showers and are larger 
also for lower primary masses. 



For Helium, the most abundant group in all QGSJet based analyses [1], the 
spectra derived from the three data sets coincide within their statistical un- 
certainties. For the three heavier mass groups no significant difference within 
systematic and statistical accuracy could be found either. 

The observed systematic differences between the solution sets are small and 
can be understood by the interplay of the dependence of shower sizes on en- 
ergy and primary particle type, increasing shower fluctuations with increasing 
zenith angle, and shifted energy threshold (caused by the fixed data range 
in A^e and A^^^) due to attenuation effects. The latter is also the reason for 
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Fig. 8. Distribution of individual Xj for the data of two different zenith angular 
ranges for the QGSJet/FLUKA solution. 

the increasing minimum primary energy of the results with increasing incli- 
nation of the incident showers. The corresponding results for the all-particle 
energy spectrum are depicted in Fig. 9. These coincide very well within their 
statistical uncertainties. 



For completeness, in Fig. 8 the Xj distributions for the two inclined shower data 
sets are shown. The Xdof parameter, which is 1.34 for vertical showers improves 
further to 1.17 for the range 18° — 25.9° and to 0.92 for the third range. This 
could be expected following the arguments above, as the QGSJet model was 
shown to have problems in describing the data at lower energies. When going 
to larger inclinations, keeping the cuts on measured particle numbers, results 
in a higher energy threshold. 
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Fig. 9. Comparison between results for the all-particle energy spectrum, using 
GHEISHA and FLUKA in the simulations. In the FLUKA case the data sets orig- 
inate from different zenith angle intervals. Error bars show the statistical uncer- 
tainty; the shaded band gives an estimate of the systematic uncertainties for the 
QGSJet/GHEISHA solution, which are of similar order for the other solutions. 
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5 Summary 



The analysis in terms of energy spectra of individual mass groups as described 
in [1] has been applied to KASCADE data based on a different low energy 
hadronic interaction model and to data sets in different intervals of the zenith 
angle. 

Using the FLUKA model instead of the GHEISHA model has no significant 
effect on the overall picture of the solution, but a slightly better description of 
the data is achieved, i.e. i.e. the correlation of electron to muon number in the 
simulations with FLUKA is in better agreement with the data which results in 
a lower Xdof distribution of the solution. The results of the unfolding analyses 
of the two-dimensional shower size spectrum for different zenith angular ranges 
show no strong or unexplainable systematic differences. Thus, the results give 
no hint to any severe problem in the simulation or the analysis, and reaffirm 
the conclusions [1] drawn from the analysis of the nearly vertical shower set: 
The knee is observed at an energy around ^ 5 PeV with a change of the index 
A7 ^ 0.4 . Considering the results of the mass group spectra, in all analyses 
an appearance of knee-like features in the spectra of the light elements is 
ascertained. In all solutions the positions of the knees in these spectra is shifted 
to higher energy with increasing element number. 

By applying the analysis to different data sets and based on different inter- 
action models, it has been demonstrated that unfolding methods are capable 
to reconstruct energy spectra of individual mass groups from air shower data, 
in addition to the all-particle spectrum. But still, the limiting factor of the 
analysis are the properties of the hadronic interaction models used and not 
the quality or the understanding of the KASCADE data. Furthermore, the 
procedure of the KASCADE data analysis, and in future also the analysis 
of KASCADE-Grande data measuring higher primary energies and muons at 
larger distances [17], gives valuable hints for the improvement of hadronic 
interaction models. The data can be confidently used when improved inter- 
action models, based on more and extended accelerator experiments, become 
available. 
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